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ABSTRACT 

We present a set of numerical simulations of stellar explosions induced by relativistic jets emanating from 
a central engine sitting at the center of compact, dying stars. We explore a wide range of durations of the 
central engine activity, two candidate stellar progenitors, and two possible values of the total energy release. 
We find that even if the jets are narrowly coUimated, their interaction with the star unbinds the stellar material, 
producing a stellar explosion. We also find that the outcome of the explosion can be very different depending 
on the duration of the engine activity. Only the longest-lasting engines result in successfiil gamma-ray bursts. 
Engines that power jets only for a short time result in relativistic supernova explosions, akin to observed engine- 
driven SNe such as SN2009bb. Engines with intermediate durations produce weak gamma-ray bursts, with 
properties similar to nearby bursts such as GRB 980425. Finally, we find that the engines with the shortest 
durations, if they exist in nature, produce stellar explosions that lack sizable amounts of relativistic ejecta and 
are therefore dynamically indistinguishable from ordinary core-collapse supernovse. 

Subject headings: Gamma-ray burst: general — Hydrodynamics — supemovae: general — supemovae: indi- 
vidual (SN2009bb) 



1. INTRODUCTION 

Massive stars end their life with powerful explosions called 
supernovae(SNe, Woosley & Weaver 1986). Some SNe are 
brighter than others, and in extreme cases a gamma-ray burst 
(GRB) - the brightest explosion in the present universe - is 
produced (Hjorth et al. 2003; Stanek et al. 2003; Woosley 
& Bloom 2006). It is believed that a GRB is associated with 
a SN when the progenitor star is compact, its stellar core has 
fast rotation, and the collapse of the core creates a fast spin- 
ning compact object that releases bipolar jets of relativistic 
matter and energy (MacFadyen & Woosley 1999; Woosley & 
Heger 2006). Recent observations (Soderberg et al. 2010; Bi- 
etenholz et al. 2010) and theoretical work (Khokhlov et al. 
1999; Burrows et al. 2007; Wheeler & Akiyama 2010; Na- 
gakura et al. 2012) have shown that an engine of the same 
kind may be present in some more mundane SNe that differ 
from their normal counterparts only by the presence of a fast 
and bright radio transient. 

Long-duration GRBs are powered by light, highly relativis- 
tic jets produced by a still mysterious compact object sitting 
at the center of an exploding massive star. The first hurdle 
that the newly born jet has to overcome to become a GRB 
is the crossing of its progenitor star. The fight, relativistic 
plasma has to bore a hole through the cold, dense stellar ma- 
terial without being excessively polluted by the star's baryons 
that would slow it down and prevent the release of the high- 
frequency photons that characterize the GRB prompt emis- 
sion (Piran 1999, Kaneko etal. 2006). Numerical simulations 
showed that the jet forms a bow shock in the stellar mate- 
rial that advances at sub-relativistic speed, shedding some of 
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its energy in the process (MacFadyen & Woosley 1999; Aloy 
et al. 2000; Zhang et al. 2003, 2004; Morsony et al. 2007; 
Mizuta & Aloy 2009). The "wasted" energy is accumulated in 
a hot cocoon that drives shock waves into the star (Ramirez- 
Ruiz et al. 2002; Lazzati & Begelman 2005; Bromberg et 
al. 2011b), unbinding it and possibly producing the super- 
nova associated with the burst (Hjorth et al. 2003; Stanek et 
al. 2003). The time it takes a typical GRB jet to cross the 
star and break out on the surface is approximately ten seconds 
(Morsony et al. 2007). However, the duration distribution 
of GRBs includes bursts as short as two seconds, requiring 
the existence of engines with a duration barely long enough 
to power the jet until the breakout time (Kouveliotou et al. 
1983). Engines with even shorter activity could exist in na- 
ture. However, their detection is challenging since they do 
not produce a clear electromagnetic signature such as a GRB. 

Observations of weak GRBs such as GRB 980425 associ- 
ated with SN1998bw (Galama et al. 1998; Kulkarni et al. 
1998), X-ray flares such as XRF080I09 (Mazzali et al. 2008; 
Soderberg et al. 2008), and relativistic radio bright super- 
novae (SN) such as SN2009bb (Soderberg et al. 2010; Bi- 
etenholz et al. 2010) call for diversity in the properties of the 
engines and/or the progenitors that power the explosion. In 
this paper we investigate the role of the duration of the en- 
gine activity on the explosive outcome of compact stars with 
a GRB-like central engine. Keeping all the jet and progenitor 
properties fairly constant, we vary the engine duration, ex- 
ploring the consequences of the non-Unearity of the jet propa- 
gation velocity for increasing engine luminosities. The paper 
is organized as follows. In Section 2 we present our numerical 
simulations, in Section 3 we describe our results, in Section 4 
we outline an analytical model to interpret our results, and in 
Section 5 we discuss the implications and limitations of our 
study. 

2. NUMERICAL SIMULATIONS 

All the simulations presented in this paper were performed 
with the FLASH code (Fryxell et al. 2000). We adopted a 
minimum resolution of 4 x 10^ cm at the highest level of re- 
finement. At this resolution the transverse dimension of the 



2 



Lazzati et al. 



injected jet is resolved into 44 elements. This resolution was 
chosen after verifying that simulations performed at twice the 
resolution would yield consistent results (Figure [Til. Our sim- 
ulations do not include magnetic fields, due to the technical 
challenge of performing MHD calculations with relativistic 
motions on an adaptive mesh. Our simulations do not con- 
sider nuclear burning since its inclusion would have made 
the simulations too long to perform at the adequate resolu- 
tion. As a consequence, ejection velocities are marginally af- 
fected. Finally, gravity from a central mass and self-gravity 
are neglected since the characteristic time-scales of the jet- 
star interaction are much shorter than the dynamical time of 
the progenitor star's collapse. Due to the absence of gravity, 
our progenitor star expands by a modest 2% of its original 
size in the 100 seconds of the simulations. The dynamical 
time scale at the stellar surface is Tdyn ~ 10^ s, while at our 
inner boundary it is Tdyn — 10 s, comparable to the break-out 
time of our jets. All our simulations adopted realistic GRB 
stellar progenitors: models 16TI and 120M from Woosley & 
Heger (2006). Model 16TI is a 16 solar-mass Wolf-Rayet star 
with an initial metallicity 1% solar and angular momentum 
L = 3.3 X 10^^ erg s. The mass of the star at pre-explosion is 
13.95 solar masses and its radius is4.1 x 10'°cm, correspond- 
ing to 0.6 solar radii. Model 120M is a 12 solar-mass Wolf- 
Rayet star with an initial metallicity 10% solar and angular 
momentum L = 2.5 x 10^' erg/s. The mass of the star at pre- 
explosion is 9.5 solar masses and its radius is 4.8 x 10'° cm, 
corresponding to 0.7 solar radii. 

In all cases, a relativistic jet with opening angle Oq = 10° 
was injected as a boundary condition at a distance = \(f cm 
from the stellar center. The engine luminosity varied among 
simulations depending on the duration of the engine activity in 
order to keep the total energy budget fixed to £ = 3 x 10^' erg 
(one set of simulations for the 16TI and 120M progeni- 
tor models) or E = 10^^ erg (one set of simulations for the 
16TI progenitor). In each simulation, the engine luminosity 
was kept constant until the cutoff time, when the engine was 
abruptly turned off. All jets were injected with Lorentz factor 
To = 5 and with enough internal energy to reach an asymp- 
totical Too = 400 upon complete, non-dissipational accelera- 
tion. This jet configuration was used numerous times to re- 
produce the properties of successful GRBs (Morsony et al. 
2007, 2010; Lazzati et al. 2009, 2011). A false color still of 
the density and expansion velocity of one of our simulations 
is shown in Figure|2] 

3. RESULTS 

The simulations were analyzed for the jet dynamics, the 
presence of an associated stellar explosion, and the nature of 
the ensuing explosion. All our simulations resulted in the ex- 
plosion of the progenitor star. However, the fraction of the 
mass that was ejected into the interstellar medium varied with 
the engine activity time. The shorter fgng was, the higher 
the fraction of the star mass that achieved escape velocity. 
Analyzing the frames at 15 s after the engine onset for our 
simulations of a the 16TI progenitor with £ = 3 x 10^' erg, 
we found that 95% of the progenitor mass had been ejected 
by our shortest-duration engine (feng = 3.0 s), while 88% of 
the progenitor mass had been ejected by our longest engine 
(feng = 15.0 s). We also found that the stellar explosions are 
not directly due to the jet, which occupies only a small frac- 
tion of the solid angle, similar to results obtained for non- 
relativistic jets (Khokhlov et al. 1999). The jet propagation. 



however, produces a hot, high-pressure cocoon that drives a 
conical shock into the stellar material. For that reason, the 
star does not explode spherically. At the time of the jet break- 
out, the star-exploding shock reaches the pole of the star but 
is still less than halfway through along the equator. For exam- 
ple, the equatorial breakout took three times longer than the 
polar one in our feng = 6.0 s simulation. Computing the frac- 
tion of ejected mass from a simulation that does not include 
gravity is non-trivial and, to some extent, not rigorous. At any 
given time in the simulation we proceeded as follows to eval- 
uate whether a parcel of matter located at radius r is bound or 
not. We first computed the amount of mass inside a sphere of 
radius r and evaluated the escape velocity from such a mass 
distribution. Subsequently, we compared the escape veloc- 
ity with the radial component of the velocity of the parcel of 
matter. If the parcel velocity is larger than the escape velocity 
we label the matter as unbound, otherwise we consider it still 
bound, eventually falling back to the stellar remnant. This 
procedure has two limitations. First the mass distribution is 
not spherically symmetric and therefore the simple approxi- 
mation that the gravity on a point mass only depends on the 
mass in the inside sphere is not correct. Second, the lack of 
gravity in the simulation affects the mass velocity itself. As 
a consequence, the values of the unbound mass fraction re- 
ported above should be considered as indicative. 

The difference in the explosive outcome is due to the non- 
linearity in the jet propagation. If the propagation velocity 
of the jet head would scale linearly with the engine luminos- 
ity, explosions with the same energy budget would look alike 
since the jet would either stall inside the star or breakout, in- 
dependently of the engine luminosity and duration. However, 
we find that the jet head propagates at a speed that scales less 
than linearly with the luminosity (Figure |3] Section 4). The 
non-linearity is brought about by the feedback between the jet 
and the high-pressure cocoon that surrounds it. A more lumi- 
nous jet has higher internal pressure that makes its head wider 
and therefore its propagation harder, shocking a larger frac- 
tion of the stellar material. This, in turn, makes the cocoon 
more energetic and the higher cocoon pressure squeezes the 
jet head facilitating its propagation. 

The different jet propagation velocities imply that in some 
configurations the jet breaks out the star's surface, while in 
others it does not. We find that there are three conditions that 
have relevance to the star's explosion outcome (Figure |4]l. If 
the engine activity time is short enough, the engine turns off 
while the jet head is still buried inside the star. As the engine 
turns off, the tail of the jet detaches from the engine and even- 
tually catches up with the jet's head. If this happens while 
the jet head is inside the star, all the bulk relativistic motion 
is lost and the star's explosion is entirely driven by the high- 
pressure cocoon that is left, resulting in an ordinary-looking 
supernova with little energy in relativistic ejecta. This condi- 
tion is shown with a dotted line in Figure]?] Alternatively, the 
engine may turn off while the jet head is mside the star but, 
by the time the tail of the jet has caught up with the head, the 
head has already broken free. In this case a small fraction of 
the material moving with bulk relativistic speed survives and, 
depending on the details, produces observable signals in the 
form of either a bright radio transient or a weak, GRB980425- 
like burst. This condition is shown in Figure [4] with a dashed 
line. Finally, if the engine active time is long enough, the jet 
breaks out while the engine is still active. In this case a fully 
developed GRB is expected. Figure |4] shows that all these 
conditions are realized by an engine with constant energy but 
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varying activity time. The comparison between three sets of 
simulations shows that the engine duration that characterizes 
each type of explosion depends also on the total energy bud- 
get and on the progenitor structure. The density and velocity 
maps of a simulation that produces the intermediate outcome 
(16 solar-mass progenitor, E = 3 x 10^' erg, t^„g = 7.5 s) are 
shown in Figure |2] 

Figure |5] shows a quantitative comparison between the re- 
sults of our simulations and the properties of type Ibc SNe 
and jet-driven stellar explosions as derived from observations 
of their associated radio transient (Soderberg et al. 2010). 
Non-thermal radio emission from SNe is produced by the in- 
teraction of the highest-velocity ejecta with the surrounding 
ambient medium (Chevalier 1998; Weiler et al. 2002; Cheva- 
lier & Fransson 2006). We selected the ejecta moving with 
(3 > 0.7 (where (3 is the velocity in unit of the speed of light) 
as the material contributing to the radio transient and we com- 
puted the expansion speed as the energy-weighted average of 
their velocity. Note that the computation was performed for 
(3 even though /3r is plotted in Figure |5] The analysis was 
performed in the last frame that we were able to compute, ei- 
ther at f = 100 s or at the time at which FLASH would run 
into a numerical instability (often occurring around the time 
of the shock breakout in the equatorial direction). The en- 
ergy in relativistic ejecta is a fairly robust number, as long as 
it is computed after the jet/cocoon breakout along the polar 
direction. Figure |6] shows the evolution of the energy in the 
ejecta of one of our longest runs as a function of time after the 
engine onset. The time of the jet/cocoon breakout along the 
polar and equatorial directions are shown. We also checked 
that the value of the energy in relativistic ejecta has a very 
small dependence on whether the actual expansion velocity 
or the asymptotic one are considered (the asymptotic expan- 
sion velocity is the one attained once all the internal energy 
has been converted in bulk outward velocity). In all cases the 
energy once the asymptotic velocity is attained is larger by 
less than 10% with respect to the one shown in the figures. 
This is consistent with the fact that the energy in Figure l6]is 
constant. Should significant acceleration take place, a marked 
evolution in the energy would be observed. In Figure [7] we 
show instead the entire cumulative distribution of the energy 
of the ejecta for three representative cases. All the simulations 
shown are performed with the 16TI progenitor and have total 
energy E = 3x 10^' erg. 

4. JET PROPAGATION MODEL 

The propagation speed of the head of the jet can be found by 
enforcing pressure balance along the discontinuity between 
the jet head and the progenitor star material. The pressure 
balance reads (Matzner 2003; Bromberg et al. 2011b): 



PjhjT'jrlifij - hf +Pj = PahaViPt+P^ 
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where /3, F, p, P, and h = I +4p/pc^ are the velocity in units 
of the speed of light, the Lorentz factor, the mass density, 
the pressure, and the dimensionless specific enthalpy, respec- 
tively. An adiabatic index 7 = 4/3 is used in the enthalpy 
expression. The subscript j refers to the jet material, the sub- 
script h refers to the jet head, and the subscript a refers to the 
ambient medium, in our case the progenitor star. The above 
equation is greatly simplified in the case of a hot, light rela- 
tivistic jet slowly advancing into a cold, high-density ambient 
medium. In that case, the pressure terms can be neglected 
and we can use the approximations F/, = 1, /3/, <c Pj, f3j = 1, 



and Pj ^ PjC^- These simplifications yield: pjhjT^j = PahaP] 
which, solving for /?/, and using the definition of h becomes: 




(2) 



Under the assumption that the pressure of the jet at the head 
is in equilibrium with the pressure of the cocoon that sur- 
rounds the jet, we have Pj = Pc = Ec/3Vc where P^, E^, and 
Vc are the cocoon pressure, energy, and volume, respectively. 
Under our assumption of sub-relativistic speed of the jet head, 
the cocoon energy is given by the energy ejected by the central 
engine: Ec = Ljt = Ljri,/c/3h, where Lj is the engine luminos- 
ity and r/, is the distance travelled by the jet head. The cocoon 
volume is given, assuming an ellipsoidal shape, by: 
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where r±^ is the transverse size of the cocoon, vj^ is the veloc- 
ity of the shock driven by the cocoon into the stellar material, 
and we have used the Kompaneets approximation in the last 
step (Begelman & Cioffi 1989; Lazzati & Begelman 2005). 
Using Eq.|3]we find the cocoon pressure as: 
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Finally, we need to find a relation between the jet pres- 
sure and Lorentz factor. The jet pressure is Pj = Ly/(4cI]yTy), 
where is the transverse cross section of the jet. In case 
of non-dissipational acceleration, the jet Lorentz factor scales 
linearly with its transverse size, and we obtain: 



4cSoP; 



1/4 



(5) 



where Eq and Fq are the initial transverse cross section and 
Lorentz factor of the jet, respectively. The system of equa- 
tions [2| |4] and[5]can now be solved for the jet head propaga- 
tion velocity: 



4r^ 
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(6) 



where we have used the identity Eq = 'nr\0\, ?"o being the jet 
injection radius and 6*0 its initial opening angle. The head 
velocity from Eq. l6]is compared to the simulation with thick 
black lines (dashea and dotted) results in Figure [3] We find 
that Eq. |6] reproduces fairly well the dependence of the head 
velocity on the jet luminosity. However, it overpredicts the 
head velocity by a factor ^ 4. This discrepancy is likely due 
to the approximation made in Eq.|5] where it was assumed that 
the jet accelerates without dissipation until it is shocked at the 
jet-star discontinuity. In fact, several recollimation shocks are 
evident in our simulations. Such shocks reduce the jet speed 
and, as a consequence, also the speed of the head (see Eq. |2]l. 

An alternative solution can be found by considering the ram 
pressure of the inner part of the jet, which is in free expan- 
sion, against the shear layer at the contact between the jet and 
star materials (Morsony et al. 2007; Bromberg et al. 2011b). 
From the results in Table 1 of Bromberg et al. (2011b) one 
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can easily find: 

The dependence of the jet head velocity on the jet luminosity 
is shown in Figure |3] with thin black lines. The results of our 
simulation are not able to distinguish between the derivation 
presented here (Eq.l6]) and the Bromberg et al. (2011b) result. 
In either case, the jetnead velocity is overestimated by a factor 
^ 4 while the dependence on the jet luminosity is qualitatively 
reproduced. 

5. DISCUSSION AND CONCLUSIONS 

We presented a set of simulations aimed at exploring the 
role of the duration of the engine activity in the explosive 
phenomenology of compact massive progenitor stars. Our 
simulations suggest that the diversity in explosive outcomes 
(GRBs, weak-GRBs, X-ray flashes, relativistic radio-bright 
SNe) is not necessarily reflected by an analogous diversity in 
progenitor and engine properties. As a matter of fact just one 
of our simulation sets, for a given progenitor structure and 
a given total energy, can give examples of the whole zoo of 
engine-driven explosions, as long as the engine activity can 
have a wide range of durations, as observed in the ORB Tgo 
distribution (Kouveliotou et al. 1983). 

Jet-driven explosions like the ones we explore, however, do 
not seem to be able to reproduce the properties of the ordi- 
nary type Ibc SNe in the lower left corner of Figure [5] This 
is expected, since stellar evolution models predict that only a 
small fraction of massive star have the properties required to 
harbor a GRB engine in their center (Woosley & Heger 2006; 
see, however, Papish & Soker 2001 for an alternative opin- 
ion). However our simulations do predict that if jet engines 
with short active times exist they produce stellar explosions. 
For example the feng = 3.0, 4.0, and 5.0 s simulations for the 
16 solar-mass progenitor and total energy E = 3 x 10^' erg 
resulted in 90% of the stellar material unbound. Such ex- 
plosions would not be associated with fast, bright radio tran- 
sients. A member of this new class of explosions may have 
been recently identified in SN2010jp (Smith et al. 2011). It 
must be noted, however, that the comparison of our simula- 
tions with observations in the lower left corner of Figure |5]is 
severely affected by the somewhat arbitrary choice of a lower 
limit /3 > 0.7 for the material to be included in the calculation 
of the "relativistic" ejecta. As shown in Figure |7] (blue line), 
the energy distribution of ordinary-looking jet-driven SNe is 
very steep and a change in the velocity threshold affects dra- 
matically the amount of energy labeled as "relativistic". In 
addition, the data from the simulations are extracted just a 
few tens of seconds after the engine onset, while the radio 
data used by Soderberg et al. (2010) were taken on a time 
scale of weeks after the burst trigger. In order to accurately 
compare the weakest radio transients in Figure [5] to our sim- 
ulations, it is required to evolve the simulations much farther 
in time and space to directly compute the radio emission, a 
task that is not feasible at the required resolution with current 
instrumentation. Additional signatures, such as asphericity, 
nucleosynthesis patterns, absorption and emission lines pro- 
files, and linear polarization need to be explored in order to 
confirm our results and find observable quantities that could 
allow us to distinguish between radio-quiet jet-driven explo- 
sions and core-collapse SNe. 

Our simulations predict that engines with a long duration 



always produce a successful burst, while only within a cer- 
tain duration interval they produce events of the intermediate 
classes, such as relativistic Ibc SNe and weak GRBs or X-ray 
flashes. Our 16TI progenitor, for example, produces inter- 
mediate events for 6 < feng < 10 s and 2 < feng < 4.5 s for 
the low-energy and high-energy engines, respectively. If the 
distribution of engine durations were flat, our mode would 
predict a larger rate of successful GRB events compared to 
events of the intermediate class. Such prediction would be 
inconsistent to current estimates, estimating the current rate 
of intermediate events either comparable to that of successful 
events (Soderberg et al. 2010) or an order of magnitude larger 
(Soderberg et al. 2006). Such comparison led Bromberg et 
al. (201 la) to conclude that low-luminosity GRBs have a dif- 
ferent origin from successful GRBs. However, the intrinsic 
engine duration distribution is not known, and it may well be 
that due to the high angular momentum that is required in the 
collapsar model to sustain prolonged accretion onto the cen- 
tral compact object (MacFadyen & Woosley 1999, Woosley 
& Heger 2006), the intrinsic distribution is skewed towards 
short durations, explaining the high rate of intermediate class 
events. Further observations and a better theoretical under- 
standing of the powering of relativistic jets will help make this 
constraint more useful for future comparison. It should also 
be pointed out that the viewing angle does have an influence 
on the burst energetics and spectrum for a given progenitor- 
engine pair (e.g., Lazzati et al. 2011), adding to the diversity 
of the events. 

Our simulations also predict a correlation between the en- 
gine duration and the energy observed as electromagnetic ra- 
diation. When trying to compare this prediction with observa- 
tions, one should keep in mind that the engine duration is not 
equal to the duration of the prompt emission (the Tgo) and that 
the energy in radiation is the "true energy", not the isotropic 
equivalent one. It is particularly difficult to connect the en- 
gine duration to the Tgo duration. In most cases of successful 
bursts one has feng > Tgo, but the difference between the two 
quantities can be large. For example, Lazzati et al. (2010) find 
that the same engine observed from different viewing angles 
can give rise to GRBs with durations spanning from a tenth 
of a second up to 100 seconds (the duration of the engine 
in that particular simulation). Such a diversity of durations 
of the radiative phase for a given engine duration is due to 
the interaction of the jet with the progenitor star. Observers 
very close to the jet axis see an initially very bright phase, 
due to the hydrodynamic collimation of the energy along the 
jet axis. Observers at larger angles, however, do not see any 
emission at early times since the jet has been collimated to 
an opening angle smaller than the angle of the line of sight. 
On the other side of the duration distribution, short engines 
can produce relatively long bursts or flashes since the dynam- 
ics of the ejecta is very different between events dominated 
by the shock breakout (e.g. Campana et al. 2006; Nakar & 
Sari 2012) and events dominated by fully developed, highly 
relativistic outflows. 

As a final remark, it should be noted that the simulations 
that we have presented here do not include all the complex 
physics of a stellar explosion. The most important limitations 
are likely the facts that magnetic fields are not considered, that 
the jet is injected at a somewhat large radius, that the simula- 
tions are performed in 2D, and that the engine/jet properties 
are not self-consistently derived from the progenitor proper- 
ties. All these limitations are likely to affect the details of 
the jet propagation. For example, a well-known 2D instabil- 
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ity produces the wedge of high-density, slow-moving mate- 
rial that is visible ahead of the jet in Figure l2] Such a wedge 
disappears in 3D simulations (Zhang et al. 2004) and its ap- 
pearance in 2D simulations has the effect of making the jet 
propagation slightly slower The bottom line is that the pre- 
cise numerical values of the engine durations that correspond 
to each class of explosions are not to be taken at face value. 
Future, more refined simulations will likely update those val- 
ues and are required to further evaluate the role of the engine 
duration in explaining the zoo of jet-driven stellar explosions. 
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Fig. 1 . — Comparison between jets from identical simulations at the adopted resolution (left panel) and at twice the adopted resolution (right panel). The upper 
part of each panel shows the density map in false colors, while the lower part shows the velocity map. The simulations shown have t^^s = 10.0 s and the frames 
are taken at f = 6.67 s after the engine onset. The two jets not only have the same size, but also show the same features, such as a big turbulent eddies on both 
sides of the head. 




Distance ( xlO^" cm) 

Fig. 2. — Density and velocity maps for the fcng = 7.5 s simulation at breakout (? = 8.13 s). The top panel shows a false-color rendering of the logarithm of the 
density, while the bottom panel shows velocity in units of the speed of light (see color scales on the right). 
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Fig. 3. — Velocity of propagation of the jet head in the stellar material at different stages of the jet-driven explosion. The red curve with squares show the 
propagation close to the core of the star, as the jet head crosses the distance r = 4 X lO' cm. The blue curve with dots show the velocity as the jet head crosses 
the distance r = 8 X lO' cm. Both datasets are obtained from the 16 solar-mass progenitor simulation with total energy £ = 3 X 10^' erg. A thick dashed line and 
a thick dotted line with equation v/, oc ]}p are overlaid on the velocity data for comparison with the analytical prediction of Eq.j6j Thin dashed and dotted lines 
show instead the comparison with the prediction of Bromberg et al. (201 1). 
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Fig. 4. — Breakout time of the jet/bubble. The black dashed Hne shows the condition for which the breakout is at the same moment at which the engine turns 
off. The black dotted line shows the condition for which the breakout takes place after the engine turns off but before the tail of the jet catches up with the jet 
head. 
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Y ■ Data from Soderberg et al. (2010) 
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H Model 16TI, E=10^^ erg, (3^ >0.7 
^ Model 120M, E=3 xio"^^ erg, >0.7 
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Fig. 5. — Comparison of the simulation results with the observational properties of Type Ibc SNe and jet-driven explosions (red symbols). The results of 
our simulations are overlaid with blue, green, and yellow symbols (circles, squares, and triangles, respectively) with the duration of the engine activity of each 
simulation indicated next to each point in the corresponding color 
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Fig. 6. — Energy in relativistic ejecta (red) and average Lorentz factor (blue) as a function of time after the engine onset for a simulation with fcng = 6.0 s and 
E = 10^' erg. The fi gure shows that measuring the energy and Lorentz factor just after the shock breakout along the polar axis yield robust results, even if the 
shock breakout along the equatorial plane has not yet taken place. 
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